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ABSTRACT 

We have estimated elemental abundances of the planetary nebula Hen2-436 in the Sagittarius (Sgr) spheroidal 
dwai-f galaxy using ESO/VLT FORS2, Magellan/MMIRS, and Spitzer/IRS spectra. We have detected candi- 
dates of fluorine [Fll]A4790, krypton [Kr lll]A6826, and phosphorus [Pll]A7875 lines and successfully esti- 
mated the abundances of these elements ([F/H]=H-1.23, [Kr/H]=H-0.26, [P/H]=H-0.26) for the first time. These 
elements are known to be synthesized by neutron capture process in the He-rich intershell during the thermally 
pulsing AGB phase. We present a relation between C, F, P, and Kr abundances among PNe and C-rich stars. 
The detections of F and Kr in Hen2-436 support the idea that F and Kr together with C are synthesized in the 
same layer and brought to the surface by the third dredge-up. We have detected N II and O II optical recombi- 
nation lines (ORLs) and derived the N^^ and O^^ abundances. The discrepancy between the abundance derived 
from the oxygen ORL and that derived from the collisionally excited line is > 1 dex. To investigate the status 
of the central star of the PN, nebula condition, and dust properties, we construct a theoretical spectral energy 
distribution (SED) model to match the observed SED with CLOUDY. By comparing the derived luminosity 
and temperature of the central star with theoretical evolutionary tracks, we conclude that the initial mass of 
the progenitor is likely to be ^1.5-2.0 Mq and the age is ~'3000 yr after the AGB phase. The observed ele- 
mental abundances of Hen2-436 can be explained by a theoretical nucleosynthesis model with a star of initial 
mass 2.25 Mq, Z=0.008 and LMC compositions. We have estimated the dust mass to be 2.9 x 10"^ Mq (amor- 
phous carbon only) or 4.0 x 10"^ Mq (amorphous carbon and PAH). Based on the assumption that most of the 
observed dust is formed during the last two thermal pulses and the dust-to-gas mass ratio is 5.58 xlO""*, the 
dust mass-loss rate and the total mass-loss rate are <3. 1 x 10"'^ Mq yr"' and <5.5 x 10"'' Mq yr"' , respectively. 
Our estimated dust mass-loss rate is comparable to a Sgr dwarf galaxy AGB star with similar metalUcity and 
luminosity. 

Subject headings: ISM; planetary nebulae; individual (Hen2-436), ISM; abundances, ISM; dust 
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1. INTRODUCTION 

Currently, >5,000 objects are regarded as planetary nebulae 
(PNe) in the local group galaxies. Of them, BoBnl (Otsuka et 
al. 2008; Otsuka et al. 2010), Wrayl6-423, StWr2-21 (Kni- 
azev et al. 2008; Zijlstra et al. 2006), and Hen2-436 (Walsh 
et al. 1997; Dudziak et al. 2000; this paper) belong to the 
Sagittarius (Sgr) dwarf galaxy. Sgr dwarf galaxy PNe are in- 
teresting objects as they provide direct insight into old, low- 
mass stars and also information on the history of the Galac- 
tic halo. In some galactic evolution scenarios, the galactic 
halo is partly built up from the tidal destruction and assimi- 
lation of dwarf galaxies. Sgr dwarf galaxy PNe are ideal lab- 
oratories and appropriate references to study the evolution of 
metal-deficient stars found in the Galactic halo and low-mass 
stars (<3.5 Mq) in the LMC. The metallicity of Sgr dwarf 
galaxy PNe is relatively low, e.g., <0.01 Zq in BoBnl or 
'--^0.5 Zq in the others, which is equal to the typical metallic- 
ity of the LMC, and all the Sgr dwarf galaxy PNe are C-rich 
([C/0]> 0; e.g., Zijlsti-a et al. 2006) based on gas-phased C 
and O abundances. The Sgr dwarf galaxy PNe and LMC PNe 
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would complement each other in a study of metal-deficient 
PNe. In addition, since the distance to the Sgr dwarf galaxy 
is well determined (24.8 kpc; Kunder & Chaboyar 2009) and 
the interstellar reddening to the Sgr dwarf galaxy PNe is rel- 
atively low {E{B-V) <0.1), we can accurately estimate in- 
trinsic flux densities. Moreover, by using the spatially highly 
resolved images taken by Hubble Space Telescope (HST) and 
8-m class ground based telescopes, we can estimate the size 
and shape of the nebulae. The nebular size is an important pa- 
rameter in building spectral energy distributions (SED). We 
have collected data on several Sgr dwarf galaxy PNe to in- 
vestigate elemental abundances, dust mass, and evolutionary 
status of the central star in a metal-deficient environment. In 
this paper, we analyze the most recently secured spectral data 
of the Sgr dwarf galaxy PN Hen2-436 (PN G004.8-22.7). 

Hen2-436 is an interesting object in terms of chemical 
abundances and dust production in metal-poor environments. 
Zijlstra et al. (2006) argued that dust exists within the neb- 
ula. Sterling et al. (2009) detected [Krlll]A2.19 fim and 
[Se IV]A2.29 /im lines in the Gemini/GNIR spectra, although 
the amounts of these elements are not estimated yet. In the 
extra-Galactic PNe, these slow neutron capture elements (s- 
process elements) have so far only been detected in BoBnl 
(Otsuka et al. 2010) and this nebula. The detection of such 
rare elements in PNe is highly interesting. 

For Hen2-436, we performed a comprehensive chemi- 
cal abundance analysis based on optical ESO/VLT FORS2 
spectra, near-IR spectra from Magellan/MMIRS, and mid-IR 
SpitzerflRS spectra. From FORS2 spectra, we found candi- 
date detections of fluorine (F), phosphorus (P), and krypton 
(Kr) forbidden lines, and we estimate abundances of these 
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elements. The estimations of F and P are done for the first 
time. Both elements are synthesized by neutron capture in 
the He-rich intershell of AGB stars. Through spectral energy 
distribution (SED) modeling with the photo-ionization code 
Cloudy (Ferland 2004), we infer the evolutionary status of 
the central star and try to estimate the dust mass. We fur- 
ther investigate the evolutionary status of the progenitors of 
Hen2-436 and the other Sgr dwarf galaxy PNe with similar 
metalUcity. The observed chemical abundances are compared 
with theoretical nucleosynthesis model predictions. 

2. DATA AND REDUCTION 
2.1. VLT/F0RS2 Archive Data 

The optical low-dispersion spectra of Hen2-436 in the range 
3320 A to 8620 A are available from the European South- 
ern Observatory (ESO) archive. The observations were per- 
formed by M.Pena (Prop. I.D.: 077.B-0430B) on August 
23th, 2006, using the visual and near UV FOcal Reducer and 
low dispersion Spectrograph 2 (FORS2; Appenzeller et al. 
1998) at the Cassegrain focus of the ANTU, one of the four 
8.2-m telescopes of the ESO Very Large Telescope (VLT) at 
Paranal, Chile. The entrance slit size was 4 17. "79 in length 
and 1"01 in width. The position angle (P.A.) was set to 0°. 
The 2x2 on-chip binning pattern was chosen, hence the sam- 
pling pitch was ~1.4 A pixel"' in wavelength and 0."5 pixel"^ 
in space. By measuring the FWHM (~4.7 A) of sky lines 
around 5500 A, we determine the spectral resolving power R 
(A/AA) to be ~1200. The data were taken by long (3 x60 sec) 
and short (3 x 1,5, and 10 sec) exposures. In order to measure 
the fluxes for strong lines such as |0 III] and Ha, we used the 
short exposure frames, and for weak line measurements we 
used the long exposure frames. The seeing was ~1."0 during 
the exposure. The standard star LTT7987 was observed for 
flux calibration and telluric absorption correction. For wave- 
length calibration, we used He-Hg-Cd lamp frames and night 
sky-lines recorded in the object frames. 

Data reduction and emission hne analysis were performed 
mainly with a long-sUt reduction package noao.twodspec in 
IRAF'. Data reduction was performed in a standard manner. 
In measuring the fluxes of emission-lines, we assumed that 
the line profiles were all Gaussian and we appUed multiple 
Gaussian fitting techniques. 

In the panels (a), (b), (c), and (d) of Figure 1 , we present the 
FORS2 spectra. The vertical axis is the scaled flux density in 
units of A"', normalized to the H/3 flux F(H/3)=1000, and the 
horizontal axis is the rest wavelength. We also present candi- 
date detections of fluorine [F ii] A4790, krypton [Kr III] A6826, 
and phosphorus [P ii]A7875 Unes. 

2.2. Magellan/MMIRS Observation 

We performed J band (1.17-1.33 fxm) moderate-resolution 

spectroscopy using the MMT and Magellan Infrared Spectro- 
graph (MMIRS; McLeod et al. 2004) at the Cassegrain focus 
of the Magellan Clay 6.5-m telescope on April 27th, 2010 
(PI.: M.Meixner). The detector of MMIRS is a 2048x2048 
pixel Hg-Cd-Te Hawaii-H array. The entrance slit size was 
420" in length and 1" in width. The P.A. was set to 0°. 
The sampling pitch is ^2.5 xlO""* /im pixel"' in wavelength 

' IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in Astron- 
omy (AURA), Inc., imder a cooperative agreement with the National Science 
Foundation. 



and 0."2 pixel"' in space. We determine R to be ^1400 by 
measuring the FWHM (^8.5 xlO""* /im) of night sky-lines 
around 1.23 jim. We observed Hen2-436 and the standard 
star HIP94663 (AOV, 2MASS 7=7.425) using a three-point 
dither pattern for sky background subtraction. The total expo- 
sure time for Hen2-436 is 900 sec (3x300 sec). The seeing 
was ^1" during the exposure. For wavelength calibration, 
we used night sky-lines recorded in the object frames. Data 
reduction was performed using IRAF. The telluric absorp- 
tion and flux caUbration were corrected using the HIP94663 
frames. The MMIRS spectrum s presented in Figure 1(e). The 
strong emission lines are He i-n Pa/3 A1.28/Ltm. 

2.3. Spitzer/IRS Archive Data 

We used the data set of program ID: P30333 (PI: A.Zijlstra) 
taken by the Spitzer space telescope on October 18th, 2006. 
The data were taken with the Infrared Spectrograph (IRS, 
Houck et al. 2004) using the SL (5.2-14.5 /im) and LL (14- 
38/im) modules. Due to the low signal-to-noise ratio, the SL 
data were not available for chemical abundance analysis and 
SED modehng. The one-dimensional spectra were extracted 
using spice version cl5.0A. We extracted the region within 
± 1 " from the center of each spectral order and summed along 
the spatial direction. No correction for interstellar extinction 
was made since it is negligibly small in this IR wavelength 
band. In Figure 1(f) we present the extracted Spitzer spec- 
trum. We detected the [Neiii] A15.5/im fine-structure line. 
The spectrum shows a broad emission feature to the red of 
the [Ne III] A15.5 /im line. Since Hen2-436 is considered 
to be a C-rich PN, the feature might be related to the C-C-C 
PAH bending modes, which are observed in some Magellanic 
Cloud PNe (Bemard-Salas et al. 2009). It should be noted, 
however, that to date, the C and O abundances of Hen2-436 
have only been derived from different types of emission lines; 
the C and O abundances are from recombination and from for- 
bidden lines, respectively. There are no UV spectra available 
from which to estimate the C^""" abundances using the C III] 
AA1906/09 Unes. Cohen & Barlow (2005) found that the PAH 
7.7 jim band is detectable in PNe with C/O > 0.56!jj;^} based 
on ISO/S WS spectra of over 40 PNe. To verify whether Hen2- 
436 has PAHs, we need to check the C/O ratio based on the C 
and O abundances from the same type of emission lines. 

2.4. HST/WFPC2 Archive Data 

The archival HST/WFPC2 images taken with the F502N 
(Ac=5013 A/AA=35.8 A), F547M (5484 A/638 A), and 
F656N (6564 A/28 A) filters are publicly available. These im- 
ages were taken by A.Zijlstra (Proposal I.D.: 9356) on May 
4fli, 2003. The WFPC2 frames reduced and calibrated by 
the HST pipeline (including MultiDrizzle), were downloaded 
from the HST archive at the Canadian Astronomy Data Cen- 
tre (CADC). In Figure 2 we present the intensity contour map 
overlaid on top of the F656N image. Hen2-436 has an elon- 
gated nebula along P.A.=90°. We performed aperture pho- 
tometry to estimate the amount of light passing through the 
slit and to do flux correction (see below). We measured the 
total flux within a 1" radius and subtracted the background 
from an annulus centered on the PN with iimer and outer 
radii of 1.1" and 1.3". We performed aperture corrections 
using point spread functions generated by the Tiny Tim soft- 
ware package^. For F547M, we assumed a 9000 K blackbody 

* http://www.stsci.edu/software/tinytim/tinytim.html 
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TABLE 1 
Observation LOGS forHen2-436 



Obs. Date Telescopes/Instruments Wavelength or Exp.Time 

Filter 



2006/08/23 


VLT/F0RS2 


3320-8620 A 


3x60,10,5,1 sec 


2010/04/27 


MageUan/MMIRS 


1.17-1.33 //m 3x300 sec 


2006/10/18 


Spitzer/IRS 


14-38 Aim 


2x14.68 sec 


2003/05/04 


HST/WFPC2 


F502N 


2x80 sec 






F547M 


1x60 sec 






F656N 


2x100 sec 



TABLE 2 

The fluxes measured from the HST images. 



Filter 



Ac 
(A) 



Flux 
(erg s~' cn 



F502N 5013 5.38(-12) ± 7.96(-14) 
F547M 5484 5.69(-13) ± 6.06(-14) 
F656N 6564 2.43(-12) ± 3.04(-14) 



function as the incident SED. The measured fluxes are listed 
in Table 2. X{—Y) means X x 10"'' hereafter. The uncertainty 
corresponds to the standard deviation of the background. To 
estimate the contributions from both the stellar and nebular 
continuum to F502N and F656N, we measured the total emis- 
sion line flux in the F547M band using the FORS2 spectra. 
Since the total emission hne flux in the F547M band is 2.34(- 
13)± 1.85(-14)ergs"i cm"^, about half of the total flux could 
originate from the continuum. Using the averaged continuum 
flux density in the F547M, we estimated the stellar and nebu- 
lar continuum subtracted [O III] and Ha fluxes to be 5.29(-12) 
and 2. 3 9 (-12) erg s"' cm"^, respectively. 

The observation logs for Hen2-436 are suimnarized in Ta- 
bid. 

3. RESULTS 

3.1. Detected Lines and Interstellar Reddening Correction 

We have detected over 100 lines, including useful lines to 
investigate plasma diagnostics and estimate ionic abundances. 
We estimate ionic abundances using detected optical recom- 
bination lines (ORLs) and coUisionally excited Unes (CELs). 

Before performing plasma diagnostics and chemical abun- 
dance analysis, we corrected our measured fluxes for inter- 
stellar reddening by determining the reddening coefficient at 
H^, c(H/3). We compared the observed hne ratio of Ha to 
H^ to the theoretical ratio computed by Storey & Hummer 
(1995) assuming a temperature T^ = 10"* K a density n^= 10^ 
cm~^ and a nebula optically thick to Ly-a (Case B of Baker & 
Menzel 1938). From this we estimate c(H/3) = 0.23 ± 0.05. 
From Seaton's (1979) relation c(H/3) = \A1E{B-V) one ob- 
tains E{B-V) = 0.16 ± 0.03, which agrees fairly well with 
the value (0. 13) measured in the direction of Hen2-436 in the 
Galactic extinction map of Schlegel et al. (1998). 

All of the line fluxes in optical to near-IR were de-reddened 
using the formula: 



log 



10 



/(A) 



am 



= log 



10 



F(X) 



Fim. 



+ c(H/3)/(A), (1) 



where /(A) and F(A) are the de-reddened and the observed 
fluxes at A, respectively, and /(A) is the interstellar extinction 
at A, computed by the reddening law of Cardelli et al. (1989) 
with Rv = 3.1. Comparing the observed [O iii]A5007, the 



Ha fluxes from the F0RS2 observation and the continuum 
subtracted F507N and F656N fluxes, we estimate 83% of the 
total flux from Hen2-436 to be passing through the 1 ."0 width 
slit entrance of FORS2. We estimate the total H/3 flux to be 
6.37(-13) ± 9.76(-15) erg s"^ cm~^. Our measured value is 
consistent wifli Cahn et al. (1992), logioF(H/3) = -12.02 ± 
0.03. Using c(H/?), we estimate the de -reddened total H/3 flux 
to be 1.08(-12) ± 1.46(-13) erg s'^ cm'l 

The de-reddened emission line fluxes detected in the 
FORS2, MMIRS, and IRS spectra are Usted in Table 3. These 
line fluxes are normalized to /(H/3)=100. In the columns of 
laboratory wavelength (Aiab), we indicate the stellar origin 
lines by the asterisks (*). These stellar lines have the wide 
FWHM, typically, > 8 A. Thanks to the large photon collect- 
ing power of VLT/FORS2, we detected 192 emission lines in 
total, including many weak intensity recombination lines of 
C, N, and O. Our measurements largely improved over previ- 
ous detections by Welsh et al. (1997), who detected 51 emis- 
sion lines in ESQ NTT 3.5-m/EMMI (spectral resolution 4.5 
A) and ESO 1.5-m/Boller and Chievens spectra (5-7 A). The 
normalized line fluxes /(A) of Welsh et al. (1997) and our 
measurements agree within 2 %. 

3.2. Plasma Diagnostics 
3.2.1. CELs Plasma Diagnostics 

We have detected several coUisionally excited lines (CEL) 
useful for estimations of the electron temperatures (r^) and 
densities (n^). We have examined the electron temperature 
and density structure within the nebula using 8 diagnostic 
CEL ratios. Electron temperature and density were derived 
from each diagnostic ratio for each hne combination by solv- 
ing for the level populations in a multi-level (> 5 for most 
ions) atomic model using the collision strengths iljj (j>i) and 
spontaneous transition probabilities Aji for each ion. The es- 
timated electron temperatures and densities are listed in Ta- 
ble 4. We used the same atomic data as in Otsuka et al. 
(2010). A diagnostic diagram that plots the loci of the ob- 
served diagnostic line ratios on the log(n^)-T^ plane is shown 
in Figure 3. The solid lines indicate diagnostics of the elec- 
tron temperature, and the broken lines are electron density 
diagnostics. This diagram shows that most of the CELs hav- 
ing ionization potentials (I.R) > 13.5 eV are emitted from 
T; ~10000-14000 K and n^ -80000-200000 cm"^ ionized 
gas. In such a high density nebula, [S ll] and [O II] diagnostic 
line ratios do not give reliable due to the low critical den- 
sity of [Sii]AA6716/31 and [On]A3727 fines. Hence, fliese 
diagnostic lines were used only for density derivation. Mean- 
while, the [N ii], [O III], [Ne III], and [Arm] diagnostic fines 
can be used as reliable indicators. 

For [N II|A5755 and |0 ll|AA7320/30, we need to consider 
the recombination contamination from N^""" and O^""" lines. 
We subtracted this contamination using Equations (2) and (3) 
given by Liu et al. (2000), 



/«([N11]A5755) 
/j;([Oll]AA7320/30) 



3.191 i^) X — 
lOV H+ 



:9.36 



104 



0.44 



i2+ 



o_ 

H+ 



(2) 
(3) 



where N^^/H"^ and O^^/H^ are the doubly ionized nitrogen and 
oxygen abundances, respectively. Adopting the value derived 
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Fig. 1 . — (panels a, h, c, and d) the FORS2 spectrum of Hen2-436. The vertical axis is the scaled flux density in unit of A which is normalized to the H/3 flux 
F(H/3)=1000. (panels e and f) the MMIRS and Spitzer spectrum of Hen2-436. The vertical axis is the scaled flux density in units of /xm"' which is normalized 
totheF(H^)=l. 

from optical recombination line (ORL) analysis (see Section 
3.3.2), we estimate /s([Nll] A5755) --0.2, which is approxi- 
mately 22 % of the observed value. Likewise, we estimated 
IrHO II] AA7320/30) >2.9, which is approximately >25 % of 
the observed value. 

First, we calculated the electron density. We adopted a con- 
stant temperature of 12 000 K for the calculations of n^(lS ii]) 
and ne([Oll]), and 14000 K when calculating nt([Cllll]) 
and ne([ArlV]). Next, we calculated the electron tem- 
perature. Tf CIN III) was calculated assuming a density of 
n,([Olll), r,(|Arlll]) was using «,([Cllll]), and r,([Nelll]) 
and T^([0 Illj) were using n^HAr IV]). The resultant and 
are Usted in Table 4. 



3.2.2. ORLs Plasma Diagnostics 

We detected optical recombination lines (ORLs) of He, C, 
N, and O. To calculate ionic abundances of these elements 
using ORLs, we estimated the electron temperature from the 
Balmer discontinuity and He I line ratios. We used the same 
atomic data employed by Otsuka et al. (2010). 

The electron temperature derived by these two methods are 



Usted in Table 4. The Balmer discontinuity electron tempera- 
ture TJB]) of 14 500 ± 2200 K was estimated using the equa- 
tion given in Liu et al. (2001). 

Using four diagnostic He I line ratios, the He I electron tem- 
perature TfiHe I) was estimated assuming a constant electron 
density = 10^ cm~^. All the He I Unes we chose for the 
TtXHe I) derivation are insensitive to electron density. We 
adopted the emissivities of He I lines from Benjamin et al. 
(1999). We consider that the r,(Hel) from the A7281/A6678 
ratio is the most rehable value because (i) He i A6678 and 
A728 1 are from levels having the same spin as the ground state 
and the Case B recombination coefficients for these lines by 
Benjamin et al. (1999) are more certain than those for A4471 
and A5876; (ii) the interstellar extinction is less of a compli- 
cation due to the similar wavelengths of these lines. Thus, we 
adopted T^iRe i) derived from the He i A7281/A6678 ratio for 
the He^ abundance calculation. 

3.3. Ionic Abundances 
3.3.1. CELs Ionic Abundances 
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TABLES 

The line list of Hen 2-436 detected in the F0RS2/MMIRS/IRS observations. 





Ion 


/(A) 


/(H/3) 




3498.77 


Hel 


0.367 


0.516 


0.371 


3587.28 


He I 


0.349 


0.659 


0.180 


3634.24 


He 1 


0.341 


0.668 


0.102 


3643.93 


Ne 11 


0.339 


0.394 


0.214 


3679.35 


HI 


0.332 


0.123 


0.110 


3682.81 


He I 


0.331 


0.182 


0.108 


3686.83 


Hi 


0.330 


0.436 


0.173 


3691.55 


Hi 


0.329 


0.643 


0.118 


3697.15 


Hi 


0.328 


0.713 


0.151 


3704.98 


He I 


0.327 


1.786 


0.134 


3711.97 


Hi 


0.325 


1.407 


0.110 


3721.94 


Hi 


0.323 


2.041 


0.158 


3726.03 


[On] 


0.322 


9.249 


0.424 


3734.37 


Hi 


0.321 


2.241 


0.157 


3750.15 


Hi 


0.317 


3.066 


0.164 


3759.88 


Oiii 


0.315 


0.351 


0.118 


3770.63 


Hi 


0.313 


3.947 


0.171 


3784.89 


He I 


0.310 


0.214 


0.147 


3797.90 


Hi 


0.307 


5.126 


0.216 


3819.60 


Hel 


0.302 


1.554 


0.451 


3835.38 


H 1 


0.299 


7.821 


0.505 


3868.77 


[Ne III] 


0.291 


53.568 


1.974 


3889.05 


Hi 


0.286 


17.481 


0.830 


3920.68 


Cii 


0.279 


0.175 


0.055 


3926.54 


Hel 


0.277 


0.181 


0.046 


3970.07 


Hi 


0.266 


30.116 


1.037 


4026.18 


Hel 


0.251 


2.564 


0.212 


4068.60 


IS II] 


0.239 


2.067 


0.245 


4076.35 


IS II] 


0.237 


0.811 


0.226 


4101.73 


HI 


0.229 


23.668 


0.726 


4120.81 


He I 


0.224 


0.479 


0.282 


4143.76 


He I 


0.217 


0.376 


0.275 


4229.27 


[Fev]? 


0.191 


0.137 


0.032 


4267.15 


Cii 


0.179 


1.102 


0.032 


4292.21 


On 


0.172 


0.040 


0.024 


4340.46 


Hi 


0.156 


46.149 


1.177 


4363.21 


[O III] 


0.149 


13.901 


0.502 


4387.93 


He I 


0.141 


0.671 


0.039 


4414.90 


On 


0.133 


0.084 


0.029 


4439.71 


[Fell] 


0.125 


0.333 


0.028 


4471.47 


He I 


0.115 


5.911 


0.150 


4571.10 


Mgi] 


0.084 


0.224 


0.021 


4634.12* 


NlII 


0.065 


0.877 


0.040 


4640.64* 


NllI 


0.063 


0.713 


0.034 


4650.25* 


C III 


0.060 


1.679 


0.050 


4658.64* 


Civ 


0.058 


1.439 


0.050 


4665.86* 


C III 


0.056 


0.948 


0.046 


4673.73 


On 


0.053 


0.522 


0.032 


4676.23 


On 


0.052 


0.339 


0.029 


4687.55 


[Fe II] 


0.049 


0.564 


0.027 


4713.17 


Hel 


0.042 


1.379 


0.028 


4740.17 


|Ar IV] 


0.034 


0.725 


0.019 


4756.59 


[Mn III] 


0.029 


0.043 


0.031 


4789.45 


[F II] 


0.020 


0.106 


0.018 


4802.58 


Ne II 


0.016 


0.088 


0.020 


4861.33 


H I 


0.000 


100.000 


2.165 


4921.93 


Hel 


-0.016 


1.649 


0.041 


4931.23 


[Gin] 


-0.019 


0.151 


0.008 


4958.91 


[GUI] 


-0.026 


265.073 


6.368 


5006.84 


[O III] 


-0.038 


792.066 


16.075 


5041.26 


On 


-0.046 


0.052 


0.033 


5047.74 


He I 


-0.048 


0.262 


0.018 


5114.06 


V 


-0.063 


0.070 


0.026 


5121.83 


C 11 


-0.065 


0.080 


0.053 



A, 



lab 



Ion 



/(A) I{HP) 5l(H/3) 



Alab 



Ion 



/(A) /(H/3) SI(HP) 



5131.25 
5145.17 
5172.34 
5179.90 
5191.82 
5197.90 
5259.06 
5261.19 
5270.40 
5335.65 
5344.31 

5376.19* 
5403.40 
5426.70 

5470.68 * 

5471.10* 
5480.05 
5488.46 
5512.77 
5517.66 
5537.89 
5543.47 
5552.68 
5577.34 
5602.44 

5592.37* 
5754.64 
5794.88 

5811.97* 
5875.62 
5889.79 
5896.78 
5910.58 
5927.82 
5931.83 
5958.38 
5977.03 
5998.71 
6046.23 
6101.79 
6118.26 
6141.70 
6151.27 
6159.12 
6300.30 
6312.10 
6346.97 
6363.78 
6461.71 
6548.04 
6562.77 
6583.46 
6678.15 
6716.44 
6730.81 
6779.94 
6791.47 
6800.68 
6826.38 
6907.80 
6933.89 
6989.46 
7001.90 
7065.18 



[Krvl? 

Cii 

N II 

C III 

[At III] 

[Nl] 

Cn 

Nn 

[Fe III] 

[Fen] 

Fein 

Cm 

Cn 

[Fe III] 

Crv 

Civ 

Nn 

Cii 

Oi 

[Cl III] 
[Cl III] 
Nil 
N II 
[O I] 
[Kvi] 
Oiii 
[N II] 
He II 
C IV 
Hel 
Cn 
Hen 
[Fen] 
Nn 
Hen 
Oi 
Hen 
[Cu III] 
Oi 
[KIV] 
Hen 
Ban? 
Cn 
[Mn III] 
|Ol| 
[S III] 
Mg II 
[O I] 
Nil 
[N II] 
Hi 
[N II] 
He I 
[S II] 
[S II] 
Cii 
Cn 
Cn 
[Kr III] 
[Fe III] 
He I 
He I 
O I 
He I 



-0.068 
-0.071 
-0.077 
-0.078 
-0.081 
-0.082 
-0.095 
-0.096 
-0.098 
-0.111 
-0.113 
-0.119 
-0.124 
-0.129 
-0.136 
-0.137 
-0.139 
-0.140 
-0.144 
-0.145 
-0.149 
-0.150 
-0.151 
-0.156 
-0.160 
-0.164 
-0.185 
-0.191 
-0.193 
-0.203 
-0.205 
-0.206 
-0.208 
-0.211 
-0.211 
-0.215 
-0.218 
-0.221 
-0.228 
-0.235 
-0.238 
-0.241 
-0.242 
-0.244 
-0.263 
-0.264 
-0.269 
-0.271 
-0.284 
-0.296 
-0.298 
-0.300 
-0.313 
-0.318 
-0.320 
-0.326 
-0.328 
-0.329 
-0.333 
-0.343 
-0.347 
-0.354 
-0.356 
-0.364 



0.081 
0.074 
0.015 
0.010 
0.063 
0.142 
0.029 
0.051 
0.079 
0.042 
0.108 
0.797 
0.054 
0.046 
0.066 
0.097 
0.046 
0.034 
0.045 
0.057 
0.160 
0.025 
0.034 
0.239 
0.074 
0.032 
0.918 
3.528 

10.124 

19.558 
0.030 
0.015 
0.008 
0.016 
0.026 
0.070 
0.034 
0.027 
0.056 
0.031 
0.030 
0.018 
0.050 
0.018 
5.523 
1.634 
0.034 
1.829 
0.159 
3.589 
285.000 

10.915 
4.923 
0.488 
1.098 
0.050 
0.067 
0.018 
0.036 
0.058 
0.016 
0.018 
0.035 

15.063 



0.022 
0.017 
0.008 
0.013 
0.005 
0.007 
0.005 
0.017 
0.011 
0.004 
0.009 
0.052 
0.006 
0.003 
0.006 
0.009 
0.003 
0.002 
0.005 
0.007 
0.014 
0.002 
0.003 
0.011 
0.009 
0.008 
0.120 
0.999 
1.300 
0.546 
0.002 
0.009 
0.005 
0.004 
0.003 
0.011 
0.009 
0.007 
0.004 
0.005 
0.006 
0.006 
0.008 
0.010 
0.197 
0.092 
0.007 
0.069 
0.010 
0.154 
11.874 
0.504 
0.194 
0.053 
0.058 
0.007 
0.011 
0.004 
0.005 
0.017 
0.007 
0.007 
0.005 
0.677 



7099.80 
7113.04 
7135.80 
7160.61 
7170.50 
7231.34 
7236.42 
7254.15 
7262.70 
7281.35 
7298.04 
7319.46 
7330.20 
7468.31 
7499.85 
7499.85 
7530.80 
7725.90* 
7751.10 
7816.14 
7816.14 
7875.99 
8046.30 
8196.50 
8203.93 
8214.50 
8236.79* 
8240.19 
8260.93 
8264.62 
8271.93 
8276.31 
8281.12 
8286.43 
8292.31 
8298.83 
8305.90 
8312.10 
8323.42 
8333.78 
8345.55 
8359.00 
8374.48 
8397.42 
8413.32 
8421.67 
8437.95 
8444.55 
8467.25 
8480.79 
8486.29 
8500.32 
8502.48 
8519.35* 
8528.97 
8545.38 
8566.92 
8581.88 
8598.39 
8616.95 
1.253([ini 
1.279Aim 
1.2S2/nTi 
15.6 fim 



[Pb u]? 
Cn 
[Ar III] 
He I 
[At IV] 
Cn 
Cn 
Gl 

[Ariv] 
He I 
He I 
[On] 
[Oil] 
Ni 
He I 
He I 
Cn 
Civ 
[Ar III] 
He I 
He I 
[Pill 
[Cl IV] 
C III 
He I 
C III 
He II 
H I 
Hi 
Hel 
Hi 
Hi 
Hi 
Hi 
Hi 
Hi 
Hi 
Cm 
Hi 
Hi 
Hi 
Hi 
He I 
He I 
H I 
He II 
H I 
He I 
H I 
He I 
He I 
|C1 nil 
H I 
He II 
Hel 
HI 
Hen 
He I 
Hi 
[Fe II] 
[Fe II] 
He I 
H I 

[Ne lll| 



-0.369 
-0.371 
-0.374 
-0.377 
-0.378 
-0.387 
-0.387 
-0.390 
-0.391 
-0.393 
-0.395 
-0.398 
-0.400 
-0.418 
-0.422 
-0.422 
-0.427 
-0.452 
-0.455 
-0.464 
-0.464 
-0.471 
-0.492 
-0.510 
-0.511 
-0.513 
-0.515 
-0.516 
-0.518 
-0.518 
-0.519 
-0.520 
-0.520 
-0.521 
-0.521 
-0.522 
-0.523 
-0.524 
-0.525 
-0.526 
-0.527 
-0.529 
-0.531 
-0.533 
-0.535 
-0.536 
-0.537 
-0.538 
-0.540 
-0.542 
-0.543 
-0.544 
-0.544 
-0.546 
-0.547 
-0.549 
-0.551 
-0.552 
-0.554 
-0.556 
-0.759 
-0.767 
-0.767 



0.075 
0.053 
7.664 
0.078 
0.035 
0.345 
0.559 
0.182 
0.043 
1.329 
0.064 
8.293 
6.746 
0.016 
0.059 
0.067 
0.101 
0.771 
1.842 
0.083 
0.083 
0.042 
0.165 
0.046 
0.062 
0.051 
0.027 
0.044 
0.019 
0.033 
0.003 
0.000 
0.019 
0.029 
0.043 
0.070 
0.103 
0.125 
0.181 
0.209 
0.249 
0.408 
0.294 
2.058 
0.328 
0.020 
0.383 
1.430 
0.428 
0.023 
0.024 
0.110 
0.428 
0.030 
0.052 
0.544 
0.005 
0.117 
0.592 
0.021 
0.442 
1.129 
16.183 
31.412 



0.006 
0.008 
0.385 
0.016 
0.006 
0.039 
0.036 
0.022 
0.021 
0.064 
0.038 
0.404 
0.354 
0.011 
0.005 
0.010 
0.012 
0.094 
0.099 
0.007 
0.007 
0.008 
0.012 
0.006 
0.011 
0.007 
0.008 
0.005 
0.004 
0.003 
0.001 
0.000 
0.009 
0.012 
0.009 
0.009 
0.011 
0.010 
0.014 
0.015 
0.017 
0.027 
0.019 
0.139 
0.058 
0.003 
0.047 
0.113 
0.028 
0.017 
0.020 
0.029 
0.037 
0.013 
0.014 
0.035 
0.006 
0.014 
0.039 
0.012 
0.133 
0.212 
1.730 
1.467 
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TABLE 4 
Plasma diagnostics. 



Parameter ID Diagnostic 


Result 




{L) [oil] (,AAD/lD/olJ/(,AA4UDy//D^ 


35 600±4900 


^cm ) 


r/~\ ttI /'\'2'70'7\//\ \'7'2'7r\/'2n\ 


71 300±3900t 




\J ) yt\l 1 V J l^A'+ / ^\J}I\A 1 i /U+A / ZD J 7 


139500±60400 




(4) [Clin](A5517)/(A5537) 


50700-149300 


Te 


(5) [N II] (AA6548/83)/(A5755) 


loeooitiioot 


(K) 


(6) [O III] (A4959+A5007)/(A4363) 


1 1 900±200 




(7) [NeIII](A15.5^m)/(A3869) 


1 1 10{)±200 




(8) [ArIII](A7135+A7751)/(A5192) 


1()2()0±400 




Hel(A7281)/(A6678) 


12 500±900 




Hei(A7281)/(A5876) 


14600±3300 




Hei(A6678)/(A4471) 


7400-9200 




Hel(A6678)/(A5876) 


11200±1200 




(Balmer Jump)/(H 11) 


14500±2200 



Corrected for recombination contribution to [O II] AA7320/30. 
' Corrected for recombination contribution to [N II] A5755. 



Fig. 2.— The intensity contour map on the HST/WFPC2 F656N (Ha) 
image. 
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Fig. 3. — Plasma diagnostic diagram. Each curve labeled 
with an ID is calculated by the following line intensity ratios; (1) 
[SlI] (AA6716/31)/(AA4069/76), (2) [Oil] (A3727)/(AA7320/30), (3) 
[ArlV] (A4740) /(A7170+A7263), (4) [ClIII] (A5517)/(A5537), (5) [Nil] 
(AA6548/83)/(A5755), (6) [OlII] (A4959+A5007)/(A4363), (7) [Nelll] 
(A15.5Atm)/(A3869), and (8) [Ami] (A7135+A7751)/(A5192). The solid 
lines indicate diagnostic lines of the electron temperature. The broken 
lines indicate diagnostic lines of the electron density. For 7'e([NlI]) and 
nj([0 II]), we corrected for recombination contributions to [N II] A5755 and 
[O ll]AA7320/30, respectively. The estimated and are summarized in 
Table 4. 

For calculating ionic abundances from the CELs, we 
adopted a three zone model for Hen2-436. Adopted and 
rie combinations for each zone and ion are Usted in Ta- 
ble 5. «e([Oll]) and ^^([Nll]) are adopted for ionic abun- 
dance derivations with LP. = 0-13.6 eV (zone 1). nt([Cllll]) 
and r<;([Arlll]) are for ions with l.P.= 13.6-28 eV (zone 2). 
ne([Ar IV]) and T^ilO III]) are for ions with LP. > 28 eV (zone 
3). 

The derived ionic abundances are listed in Table 6. In the 
last line of the line series of each ion, we present the adopted 
ionic abundances in bold face. These values are obtained by 
the hne intensity weighted mean when we detected two or 
more lines. We estimated 14 ionic abundances by solving a 
>5 level atomic model. 

We detected candidates of [Fii]A4789, [Pii]A7875, and 
[Kriii]A6826 as earUer described. The estimations of the 



TABLE 5 

Adopting T, and u, i-or CELs ionic abundanci; calculations. 



Zone Ions 



ne(cm ^) 



r,(K) 



1 N+, 0+, P+, S+ Me([Oii]) 7;([NII]) 

2 F+, S-+, Cl2+, Ar2+, nj([Cllll]) ^^([Arni]) 

3 0^+, Ne^+, Cl^+, Ar3+ >ie([Arlv]) r;([Olll]) 

F^, P^, Fe^^, and Kr^'^ abundances are probably done for the 
first time. For [P ll], we adopted the transition probabili- 
ties of Mendoza & Zeippen (1982), the collisional impacts of 
Tayal (2004), and the level energy listed in Atomic Line List 
v2.05bl2^ For [Kr III], we adopted the transition probabili- 
ties of Bi6mont & Hansen(1986), the collisional impacts of 
Schoning (1997), and the level energy listed in Atomic Line 
Listv2.05bl2. 

The enhancements of these elements would give constraints 
to parameters (the third dredge-up efficiency, '^C pocket 
mass, and the number of thermal pulse, etc.) in nucleosynthe- 
sis models of low-mass stars. The detection of these elements 
is therefore highly interesting. 

For the and estimations, we subtracted the recom- 
bination contamination to [Nll]A5755 and [O ll]AA7320/30 
from N'^"'" and O^^ lines, respectively. Hence, we obtained 
the appropriate abundances, i.e., the amounts derived from 
their nebular lines. S"^ abundances from [S ii] A4068/76 lines 
are lower than those from [S ll] A67 17/31 in adopted electron 
density. S"^ abundances derived from [S ii] A4068/76 fines 
give more reliable values than those from [S ii] A67 17/31. To 
estimate Fe^^ abundance, we solved a 33 level model (from 
5Z)3tob3/'2)for[Feiii]. 

3.3.2. ORLs Ionic Abundances 

The estimated ionic abundances derived from ORLs are 
listed in Table 7. This work may constitute the first ever es- 
timation of the N and O ORL abundance. The ORL ionic 
abundances are derived from 



A^(X"+) _ a{HB) A(X'"+) /(X™+) 
NiU+) ~ a(X'"+) A(H/3) /(H/3) ' 



(4) 



' see http://www.pa.uky.edu/ peter/newpage/ 
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TABLE 6 
CEL IONIC ABUNDANCES. 







Transition 


/(A) 


SI(X) 




(5X'"+/H+ 




(A///m) 


(lower - 


- upper) 










N+ 


5754.6 






7.17(-1) 


1.35(-1) 


3.49(-6) 


1.38(-6) 




6548.0 




-•Z52 


3.59(0) 


1.54(-1) 


3.45(-6) 


7.91 (-7) 




6583.5 


3ft- 


-'02 


1.09(-i-l) 


5.04(-l) 


3.54(-6) 
3.52(-6) 


8.16(-7) 
8.37(-7) 


0+ 


3726.0/28.8 


' ''S3/2-^-D3/2,5/2 


9.25(0) 


4.24(-l) 


3.02(-5) 


1.01 (-5) 




7319.0 




■A/2,3/2 


6.39(0) 


3.11(-1) 


3.58(-5) 


1.50(-5) 




7330.0 


"'^3/2-^/2,3/2 


5.20(0) 


2.73(-l) 


3.54(-5) 


1.48(-5) 














3.32(-5) 


1.28(-5) 


q2+ 


4363.2 


'Dr 




1.39(+1) 


5.02(-l) 


1 90(-4) 


1.83(-5) 




4931.8 




-ID2 


1.51(-1) 


8.35(-3) 


2.71(-4) 


2.01(-5) 




4958.9 






2.65(-i-2) 


6.37(0) 


1.86(^) 


1.01(-5) 




5006.8 


3P2- 




7.92(+2) 


1.61(+1) 


1.93(^) 
1.91(-4) 


1.02(-5) 
1.03(-5) 




4789.5 


3P2- 


-'02 


1.06(-1) 


1.82(-2) 


7.19(-8) 


1.50(-8) 


Ne2+ 


3868.8 




.1D2 


5.36(-i-l) 


1.97(0) 


3.28(-5) 


2.22(-6) 




15.6 


3P2- 


-3P, 


3.14(+1) 


1.47(0) 


3.98(-5) 
3.54(-5) 


1.92(-6) 
2.10(-6) 


P+ 


7875.0 


'D2- 


i^o 


4.21(-2) 


7.69(-3) 


2.31(-8) 


7.37(-9) 


S+ 


4068.6 


4c. 
■J3/2" 


- P3/2 


2.07(0) 


2.45(-l) 


2.22(-7) 


6.17(-8) 




4076.4 


'Syr 


M/2 


8.11(-1) 


2.26(-l) 


2.63(-7) 


9.89(-8) 




6716.4 


4c 
■^3/2- 




4.88(-l) 


5.35(-2) 


4.02(-7) 


9.16(-8) 




6730.8 


''S3/2- 


-'D3/2 


1.10(0) 


5.82(-2) 


4.12(-7) 
2.96(-7) 


8.59(-8) 
7.77(-8) 


S2+ 


6312.1 


'O2 


-'So 


1.63(0) 


9.17(-2) 


4.59(-6) 


7.11(-7) 


CP 


5517.7 


■53/2- 


-'-D5/2 


5.70(-2) 


6.81(-3) 


5.33(-8) 


8.42(-9) 




5538.0 


^52-3/2^3 


1.60(-1) 


1.41(-2) 


5.40(-8) 


7.37(-9) 














S.38{-8) 


7.65(-9) 


a3+ 


8046.3 


3P2- 


-'Z)2 


1.65(-1) 


1.16(-2) 


8.11{-9) 6.30{-10) 




5191.8 


'O2 


-'So 


6.34(-2) 


5.46(-3) 


6.78(-7) 


1.30(-7) 




7135.8 


3^2- 




7.66(0) 


3.85(-l) 


6.71 (-7) 


6.52(-8) 




7751.1 


3Pl- 


->D2 


1.84(0) 


9.90(-2) 


6.72(-7) 
6.71(-7) 


6.66(-8) 
6.59(-8) 




4740.2 


'S3/2- 


-^i)3/2 


7.25(-l) 


1.89(-2) 


8.58(-8) 


4.22(-9) 




7170.5 


'O3/2- 


-'/'3/2 


3.47(-2) 


6.48(-3) 


8.71(-8) 


1.72(-8) 




7262.7 


'£'3/2- 


-'/'1/2 


4.29(-2) 


2.10(-2) 


1.24(-7) 
8.79(-8) 


6.10(-8) 
7.81(-9) 


Fe2+ 


5270.4 


5D3- 


3P22 


7.88(-2) 


1.08(-2) 


7.91(-8) 


1.36(-8) 




6826.4 


3P2- 


->D2 


3.57(-2) 


4.83(-3) 


3.11(-9) 5.00(-10) 



Note. — We corrected for recombination contributions to [N 11] A5755 
and [On] AA7320/30. 



where a(X™'^) is the recombination coefficient for the ion 
X'"^. From the logn^-i; plot of the [On] A3727/AA7320/30 
ratio and estimated values of T^iBJ) and T^iHei), we calcu- 
late the electron density to be ^10^ cm~^. We adopted this 
value for all ORL ionic abundance calculations. We adopted 
T^(BJ) for all recombination fines except for He"'". We adopted 
T^Qie I) for He"'" abundances. 

Effective recombination coefficients for the lines' parent 
multiplets were taken from the references listed in Table 11 of 
Otsuka et al. (2010). The recombination coefficient of each 
line was obtained by a branching ratio, B(A,), which is the ra- 
tio of the recombination coefficient of the target line, a(A;) 
to the total recombination coefficient, J2i Cfi^t) ^ ^ multiplet 
line. To calculate the branching ratio, we referred to Wiese 
et al. (1996) except for the Oil 3d-4-f transition line Oil 
A4292.2. For this line, the branching ratios were provided by 
Liu et al. (1995) based on an intermediate coupUng. We ap- 
pUed the Case B assumption (Baker & Menzel 1938) for lines 
of levels with the same spin as the ground state and the Case A 
(optically thin in Ly-a) assumption for the other multiplicity 
fines. In the last line of the line series of each ion, we present 



TABLE 7 
ORL IONIC ABUNDANCES. 



X'""^ Aiab /(A) 5/(A) X^W (5X'«+/H+ 

(A) 

He+ 6678.2 4.29(0) 1.69(-1) 1.06(-1) 1.06(-2) 

4471.5 4.91(0) 1.24(-1) 1.33(^1) 1.36(-2) 
4387.9 6.24(-l) 3.58(-2) 1.06(-1) 1.31(-2) 
4921.9 1.50(0) 3.71(-2) 1.09(-1) 9.93(-3) 

5875.6 1.42(+1) 3.96(-l) I.12(-l) 1.29(-2) 

1.15(-1) 1.25(-2) 

4267.2 1.10(0) 3.16(-2) 1.14(-3) 2.65(^) 

6151.3 5.05(-2) 8.44(-3) 1.14(-3) 2.99(-A) 

7231.3 3.45(-l) 3.86(-2) 8.04(-4) 2.47(-4) 

7236.4 5.59(-l) 3.56(-2) 6.52(-^) 1.91(-^) 

9.50(-^) 2.43(-^) 

C3+ 8196.5 4.57(-2) 5.60(-3) 6.97(-5) 1.72(-5) 

5480.1 4.56(-2) 3.24(-3) 5.69(-4) 1.22(^) 
5927.8 1.57(-2) 3.81(-3) 5.30(^) 2.20(^) 

5.59(-4) 1.47(-4) 

0^+ 4292.2 3.98(-2) 2.36(-2) 1.61(-3) 1.01(-3) 

4676.2 3.39(-l) 2.94(-2) 3.30(-3) 6.81(^) 
3.12(-3) 7.15(-^) 

Note. — For He"*", the contribution of coUisional excitation was sub- 
tracted. Oil A4291.26,4291.86,4292.21,4292.98 are observed as a blended 
line (See the text). 

the adopted ionic abundance and the error estimated by the 
weighted mean of the line intensity. In the following, we give 
short comments to several ionic abundance estimations. 

The He""" abundances were estimated using 5 different He I 
lines selected to be insensitive to electron density in order to 
reduce intensity enhancement by collisional excitation from 
the He*^ 2s level. The collisional excitation from the He'' 
2s ^5 level mainly enhances the intensity of the triplet Hei 
lines. We removed this contribution using the formulae given 
by Kingdon & Ferland (1995). The collisional excitation con- 
tamination for our measured lines was estimated to be 14.8 % 
for 6678 A, 20.5 % for 4471 A, 7.5 % for 4388 A, 9.7 % for 
4921 A, and 37.9 % for 5876 A. The derived He"" abundances 
well agreed with each other within error. 

We detected several C III lines. Of them, we chose C III 
A8196 to estimate the C^"'" abundance. Its line width is narrow 
(3.97 A), therefore this line could be of nebular origin. Since 
the ground state of C iii is singlet (2s^ ^S), we adopted the 
Case A assumption for this line. We detected several C IV and 
O III lines of wide line width. Since they would be of stel- 
lar wind origin, we did not estimate their ionic abundances; 
Hen2-436 could not be a highly ionized nebula, because the 
high I.P. nebular fines such as [Neiv] were unseen in the 
FORS2 spectra. The central star of the PN is not hot enough 
to ionize these lines. 

We measured O II A4292.2 (doublet, 3<i-4/) and A4676.2 
(quadruplet, 3s-3p) lines as well. 11 A429 1 .2 is blended with 
On AA4291.86,4292.98 A. Since the ground level of On is 
quadruplet, we adopted Case A for the doublet and Case B 
for the quadruplet lines. We found that the discrepancy be- 
tween O^"^ ORL and CEL abundances is >1 dex. The O^^ 
discrepancy has been found in many PNe including BoBnl 
(see section 4 of Otsuka et al. 2010). 

3.4. Elemental Abundances 

To estimate elemental abundances, we estimated the unob- 
served ionic abundances using an ionization correction factor, 
ICF(X). lCFs(X) for each element are listed in Table 8. 

The He abundance is the sum of He"'" and He*'^ abundances. 
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and we estimated the unseen He abundance. The C abun- 
dance is the sum of the C"^, C^"^, C^^ abundances. We cor- 
rected for the unseen C"^ assuming (C^/C) = (N^/N)cels- The 
N abundance is the sum of and N^"^. For the CEL N abun- 
dance, N^"^ was unobserved and calculated using Table 8. For 
the ORL N abundance, we accounted for the unseen N""" as- 
suming (N/N^)oRLs = (Ar/Ar^^)cELs- The O abundance is the 
sum of the O"*" and O^""" abundances. For the ORL O abun- 
dance, we assumed (0^'^/0)orls = (0^'^/0)cels- The Ne abun- 
dance is the sum of the Ne"^ and Ne^"^ abundances, and we 
corrected for the unseen Ne"^. Assuming that the F abundance 
is the sum of the F^ and F^""" abundances, we estimated the un- 
observed F^^ abundance. The S abundance is the sum of the 
S""", S^""", and S^"*" abundances. We accounted for the unseen 
abundance using the CEL O and O"^ abundances. We assume 
that the CI abundance is the sum of Cl^, Cl^"^, and Cl'"^ abun- 
dances. The unseen CI"*" is accounted for by assuming (Cl/CV') 
= (0/0"'")cELs- The Ar abundance is the sum of the Ar^""" and 
Ar^"^ abundances, and we accounted for the unseen Ar"^. The 
Fe abundance is the sum of the Fe^"^ and Fe^"^ abundances, and 
we accounted for the unseen Fe^"^. The Kr abundance is the 
sum of the Kr"*", Kr^""" and Kr^""" abundances, and we accounted 
for the unseen Kr""" and Kr^""". 

The resultant elemental abundances are Usted in Table 9. 
In the third colunon, the number densities relative to H are 
listed. In the fourth column, the logarithmic number densi- 
ties are given when log[o(H) is 12. In the fifth, sixth, and 
seventh columns, we present the enhancements relative to 
the solar abundances. For the solar abundances, we adopted 
the amounts recommended by Lodders (2003). In the last 
two columns, we present the number density in the form of 
logio(X/H)H-12 estimated by Walsh et al. (1997) and Dudziak 
et al. (2000). The amounts of Walsh et al. (1997) were es- 
timated based on optical spectra and those of Dudziak et al. 
(2000) by a photo-ionization model. Our estimated amounts 
agree fairly well with their values except N and S. The over- 
abundances of these elements relative to Walsh et al. (1997) 
would be due to the adopted for N""", S"^ and S^"'"estimations. 
Walsh et al. (1997) adopted 12600 K for all ionic abundance 
estimations, while we adopted 10600 K for N"^ and S"^. If 
we adopted a of 12 600 K, our N and S abundances would 
decrease by -0.17 and -0.54 dex, respectively. 

Our estimated [Cl,Ar/H] are almost consistent within the 
estimated error. Since both elements are not synthesized in 
low-mass stars and do not couple with the dust grains, they 
would indicate the metalUcity of the progenitor. The metal- 
licity of Hen2-436 (^^-0.5) is comparable to a typical LMC 
metallicity. The large depletion of the Fe abundance ([Fe/H]) 
indicates that most of the Fe exists not as ionized gas but as 
dust grains. 

The discrepancy between the ORL and CEL O abundances 
(>1 dex) is not due to the adopted electron temperatures; the 
adopted for the ORL O^^ estimations is consistent with that 
for CEL O^""". The O ORLs might be emitted from 0-rich and 
H-deficient knots, similar to those observed in Abell 30 (e.g.. 
Wesson et al. 2003) because the central star of Hen2-436 is a 
[WC] star and most stars of this type are H-deficient. Accord- 
ing to the classification by Acker & Neiner (2003), the central 
star of Hen2-436 is [WC4-6] based on the instrumental broad- 
ening subtracted FWHM of C iv A5812 (29.6 A) and the Une 
ratios of C ii,iii to C iv A5812. Girard et al. (2007) classified 
Hen2-436 as [WC4]. In this paper, we will not examine the O 
discrepancy problem further. To resolve the O abundance dis- 
crepancy, it is necessary to obtain high-dispersion spectra to 



increase the chance of O II detections and properly de-blend 
these O II lines with the others. 

We are able to estimate the C/O ratio of Hen2-436 to be 
0.54 ± 0.41 from the elemental ORL C and O abundances. 
This is the first time the C/O ratio has been calculated for this 
nebula using the same type of emission line. In Hen2-436, the 
CEL C abundance has not previously been estimated as men- 
tioned earlier. So far, both the CEL and ORL C abundances 
in PNe have been estimated by Wang & Liu (2007), Wesson 
et al. (2005), Liu et al. (2004), and Tsamis et al. (2004). 
From these works, we found the C2+(ORL)/C2+(CEL) of 4.10 
± 0.49 among 58 PNe and the C^+(OLR)/C^+(CEL) of 4.52 
± 0.85 among 14 PNe. Using those factors, the CEL C^+ and 
C^"*" abundances are extrapolated to be 2.32(^) ± 6.34(-5) 
and 1.54(-5) ± 4.78(-6), respectively. The expected CEL C 
abundance is 2.89(-4) ± 1 .95(-4) and the CEL C/O is 1 .68 ± 
1.08 when we adopt the same ICF listed in Table 8. 

The F, P, and Kr abundances were estimated. These ele- 
ments are known to be synthesized by neutron capture pro- 
cesses in the He-rich intershell region during the thermally 
pulsing AGB phase (e.g., Abia et al. 2010 and Lugaro et al. 
2004 for F production, Herwig 2005 and Busso et al. 1999 
for AGB nucleosynthesis including the s-process). Concern- 
ing low-mass stars (<4 Mq), most of the free neutrons are 
released by the '^C(a,n)'^0 reaction in the interpulse phase. 
^^C isotope is produced by partial mixing of the bottom of the 
H-rich convective envelope into the outermost region of the 
'^C-rich intershell layer. The nucleosynthesis models for low- 
to intermediate mass stars predict that C and neutron capture 
elements produced in the He-rich intershell are brought up to 
the stellar surface by the third dredge up. Figure 4 (a), (b), and 
(c) shows the relations between [C/Ar] and [F/Ar], [P/Ar], 
and [Kr/Ar] among PNe and C-rich stars, respectively. The 
trends between [C/Ar] and [F,P,Kr/Ar] support the theoretical 
description of the i-process and the third dredge-up occurring 
in AGB stars. 

Among Sgr dwarf galaxy PNe, F lines have been detected 
in the halo PN BoBnl (Otsuka et al. 2010, 2008), which is 
also suspected to be an old Sgr dwarf galaxy PN. Hen2-436 
would be the second detection case. 

The only stable P isotope is ^^P. ^^P is synthesized in the 
He-rich intershell via the following reaction. 

2«Mg(n, 7)2^Mg(^+)2^Al(n, 7)'^A1(^+) 

^^Si(n,7)2'Si(n,7)^°Si(n,7)3^Si(^"')^ip (5) 

In Table 10, we list the P abundances estimated for 9 
PNe. The listed P abundances are of nebular origin and 
they are estimated using [P ll]AA4669,7875 (NGC6572, 
NGC6741, NGC7027, and Hen2-436), [Pii]AA1.14/1.19 
urn (IC5117 and NGC7027), [Piii]A17.89 /zm (NGC40, 
NGC2392, NGC6210, and NGC6826). We did not detect 
[P II] A A 1 . 1 4/1 . 1 9 ^m in the MMIRS spectrum due to the low 
S/N. 

The 2.5 Mq models of Karakas et al. (2009) predict that the 
P abundance is not largely enhanced after AGB nucleosynthe- 
sis (<0.1 dex) with or without a partial mixing zone. The P 
production in the He-rich intershell depends on the amount 
of extra mixing included in the calculations to produce the 
'^C pocket (Werner & Herwig 2006). According to Werner & 
Herwig (2006), P has been detected in a number of PG 1159 
stars by the identification of the P V AAl 118/28 lines, and 
their P abundances are roughly solar. PG 1159 stars are H- 
deficient. Marcolino et al. (2007) observed 4 PNe with WR 
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type central stars using the far-UV spectrograph FU SE, and 
they found P V lines in BD+30°3639, NGC40, and NGC5315. 
They argued that the P abundance in NGC5315 is 4-25 times 
the solar abundance and the value agrees with nucleosynthesis 
calculations of Werner & Herwig (2006). For the other PNe, 
they did not give the P abundances. Our estimated [P/H] abun- 
dance could agree with Karakas et al. (2009) and Werner & 
Herwig (2006). However, except for Hen2-436 and NGC40, 
PNe listed in Table 10 do not have WR central stars, and 4 
PNe of these no-WR 7 PNe show >0.2 dex P enhancement. 
To verify the discrepancy of P enhancement between the ob- 
servation and the models, we need to increase the detection 
cases, for WR-type PNe in particular. We should note that 
the production of P in AGB models is strongly dependent on 
a number of modeUng uncertainties including reaction rates, 
the adopted convective mixing model, and the mass-loss rate; 
for example, Karakas & Lattanzio (2007) and Karakas (2010) 
adopted different reaction rates from each other, and the mod- 
els by the former produced about a factor of 10 more P than 
those of the latter. The P abundance could give a constraint 
to uncertain parameters in the AGB nucleosynthesis models, 
including the '■'C pocket mass, and it would be an efficient 
diagnostic tool, because P lines can be observed from the far- 
UV to the mid-lR. 

Among PNe in which F, P, and Kr have been detected, 
Hen2-436 seems to be very similar to NGC40 in elemen- 
tal abundances, dust composition, and central star properties. 
The elemental abundances of NGC40 are hsted in the last 
column of Table 9. In terms of elemental abundances, the 
C/O ratio and the [Cl,Ar,Fe/H] abundances are comparable to 
those of Hen2-436 as present in Table 9. The ORL and CEL 
C/O ratios of NGC40 are 0.46 and 1.41, respectively The 
C(ORL)/0(CEL) ratio is 8.32 (7.25 ± 4.86 in Hen2-436). F, 
P, and Kr abundances are also comparable. In addition, the en- 
hancements of ORL N and O abundances are relatively simi- 
lar to those of CEL abundances. Liu et al. (2004) reported that 
the C2+(ORL)/C2+(CEL) and 02+(ORL)/02+(CEL) ratios are 
5.8 and 17.8 (16.3 ± 3.9 in Hen2-436), respectively. Ramos- 
Larios et al. (2010) found PAHs emissions in the ISO/SWS 
and SpitzerflRS spectra of NGC40. Since our estimated CEL 
and ORL C/O ratios of Hen2-436 are larger than the PAH 7.7 
/zm band detection limit by Cohen & Barlow (2005) and the 
chemistry of Hen2-436 seems to be similar to NGC40, Hen2- 
436 might have PAHs. The central stars of both PNe have 
wide FWHM C iii, C iv, and N iii stellar origin lines, and they 
have WR-type central stars. De Marco & Barlow (2001) clas- 
sified NGC40 as [WC8]. The strong wind could contribute to 
enhancements of C, F, P, and Kr abundances in the nebula. 

4. DISCUSSION 

4.1. Photo-Ionization Modeling 

To investigate the properties of the ionized gas, dust, and the 
central star of the PN in a self-consistent way, we constructed 
a theoretical photoionization (P-1) model which matches the 
observed flux of emission fines and the spectral energy distri- 
bution (SED) between optical and nud-IR wavelengths, using 
Cloudy c08.00 (Ferland 2004). 

To construct an accurate P-1 model, we require information 
about the distance from us to the object, the incident SED 
from the central star and the elemental abundances, geometry, 
density distribution, and size of the nebula. Since Hen2-436 
is in the core of the Sagittarius dwarf galaxy as we mentioned 
in the introduction, we fix the distance at 24.8 kpc (Kunder 
& Chaboyer 2009). Dudziak et al. (2000) estimated the ef- 



TABLE 8 

Adopted Ionization Correction Factors (ICFs). 



X 


Line 


ICF(X) 


X/H 


He 


ORLs 


(S++S^+)/S^+ 


ICF(He)He+ 


C 


ORLs 


(1-(N+/N)r' 


ICF(C)(C2++C^+) 


N 


CELs 


(0/0+) 


ICF(N)N+ 




ORLs 


(Ar/Ar2+) 


ICF(N)N2+ 





CELs 


1 


0++02+ 




ORLs 


(0/02+)cELs 


ICF(0)02+ 


F 


CELs 


(0/0+) 


ICF(F)F+ 


Ne 


CELs 


(0/02+)cELs 


ICF(Ne)Ne2+ 


P 


CELs 


(S/S+) 


ICF(P)P+ 


S 


CELs 


[l-(I-(0+/0))3]-'/' 


ICF(S)(S++S2+) 


CI 


CELs 


(l-(0+/0))-i 


ICF(Cl)(Cl2++Cl3+) 


At 


CELs 


(l-(N+/N))-i 


ICF(Ar)(Ar2++Ar3+) 


Fe 


CELs 


(0/0+) 


ICF(Fe)Fe2+ 


Kr 


CELs 


C1/C12+ 


ICF(Kr)Kr2+ 



fective temperature T^a = (7.0±1.0)x 10"* K and luminosity 
L*=(5.4±0.4)x 10^ Lq of the central star in their P-1 model. 
Guided by their Teff and L*, we used a series of theoretical 
atmosphere models with a range of values of Jeff to supply 
the SED from the central star. We used Thomas Rauch's non- 
LTE theoretical atmosphere model'" for halo stars ([X,F]=0, 
[ZJ=-1) and we considered the surface gravity log ^=6.0 and 
6.5 cases. We varied Ten and to match the observations. For 
the elemental abundances NiXyNiH), we used the observed 
values fisted in Table 9 as a first guess. For initial N and O 
abundances, we adopted the values derived from the CELs. 
Note that line fluxes of ORLs from N and O calculated by 
the models would be underestimated. We did not consider the 
F, P, and Kr abundances because Cloudy does not calculate 
[Fii]A4790, [Kriii]A6826, [Pii]A7875 and also because very 
small abundances of these elements are not effective in cool- 
ing. For F, P, and Kr and the other unobserved elements, we 
therefore fixed [X/H]=-0.6. We fixed the outer nebular shell 
^out at 0.26", which corresponds to the FWHM of the HST 
image (Figure 2). We adopted a hydrogen density (A^h) 
profile. We also varied N{X)/N{}i), Nu at the iimer radius, 
and the inner radius Rm over a small range of values to match 
the observed line fluxes detected in the F0RS2, MMIRS, and 
IRS and 2MASS JHKs bands and our interesting mid-infrared 
band. IRS B is the integrated flux between 17 and 23 fxm and 
IRS C is between 27 and 33 /xm. 

Dust grains co-exist along with the gas in the nebula of 
Hen2-436. As discussed in Sections 2.3 and 3.4, Hen2-436 
might have PAHs. Therefore, we considered two dust com- 
position models; i) amorphous carbon (am C, hereafter) only 
(model A hereafter) and ii) am C + PAH grains (model B, 
hereafter). The optical constants were taken from Rouleau & 
Martin (1991) for amorphous carbon and taken from Desert 
et al. (1990), Schutte et al. (1993), Geballe (1989), and Breg- 
man et al. (1989) for PAHs. In the Cloudy model we as- 
sumed that the gas and dust co-exist in the same sized ionized 
nebula, so that we consider the warm dust (> 100 K) only. We 
adopted a standard MRN a~^^ distribution (Mathis, Rumpl & 

see http://astro.uni-tuebmgen.de/~rauch/ 
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Fig. 4. — (left panel): The [F/Ar]-[C/Ar] diagram. The typical error is indicated by the cross (the same for the other diagrams) . The data of BoBnl are from 
Otsuka et al. (2010). The references for Galactic PNe and MS/S stars are in Otsuka et al. (2008). {middle panel): The [P/Ar]-[C/Ar] diagram. The references for 
Galactic PNe are in Table 10. {right panel): The [Kr/Ar]-[C/Ar] diagram. The data for BoBnl (upper limit) is from Otsuka et al. (2010). The data for Galactic 
PNe are taken from Sterling & Dinerstein (2008). 



TABLE 9 

Elemental abundances of Hen2-436. 



X 


Line 


X/H 


logfX/H)+12 


fX/Hl 


[X/Arl 


fX/Olt 


Rcf.d) 


Rcf.(2) 


Ref.(3) 








(H=I2) 








Hcn2-436 


Hcn2-436 


NGC40 


He 


ORLs 


1.22(-l)±1.35(-2) 


11.09±0.05 


+0.19±0.05 


+0.79±0.33 


+0.53±0.08 


11.02±0.02 


11.03±0.01 


11.08 


C 


CELs 
















8.84 




ORLs 


1.20(-3)±7.97(^) 


9.08±0.35 


+0.69±0.35 


+ 1.29±0.48 


+1.03±0.36 




9.06±0.09 


9.61 


N 


CELs 


2.37(-5)±1.09(-5) 


7.38±0.22 


-0.45±0.24 


+0.1 5 ±0.41 


-0.11 ±0.25 


6.97±0.16 


7.42±0.06 


7.93 




ORLs 


7.42(^)±5.10(^) 


8.87±0.37 


+1.04±0.38 


+1.64±0.51 


+1.38±0.39 






9.32 





CELs 


2.24(^)±1.64(-5) 


8.35±0.03 


-0.34±0.06 


+0.26±0.34 


+{).{)0±0.08 


8.29±0.08 


8.36±0.06 


8.69 




ORLs 


3.67(-3)±9.04(^) 


9.56±0.11 


+0.87±0.12 


+1.47±0.35 


+ 1.21±0.13 






9.95 


F 


CELs 


4.85(-7)±2.15(-7) 


5.69±0.21 


+1.23±0.22 


+1.83±0.39 


+1.57±0.22 






5.48 


Ne 


CELs 


4.15(-5)±4.51(-6) 


7.62±0.05 


-0.25±0.11 


+0.35±0.35 


+0.09±0.13 


7.57±0.08 


7.54±0.06 


8.01 


P 


CELs 


5.26(-7)±2.38(-7) 


5.72±0.21 


+0.26±0.22 


+0.86±0.39 


+0.60±0.22 






5.38 


S 


CELs 


6.73(-6)±1.23(-6) 


6.83±0.08 


-0.36±0.09 


+0.24±0.34 


-0.02±0.11 


6.30±0.08 


6.59±0.05 


6.41 


CI 


CELs 


7.26(-8)±1.48(-8) 


4.86±0.09 


-0.47±0.11 


+0.13±0.35 


-0.13±0.12 






4.91 


Ar 


CELs 


8.91(-7)±5.59(-7) 


5.95±0.32 


-0.60±0.33 


+0.00±0.47 


-0.26±0.34 


5.76±0.12 


5.78±0.08 


5.94 


Fe 


CELs 


5.34(-7)±2.28(-7) 


5.73±0.20 


-1.74±0.20 


-1.14±0.39 


-1.40±0.21 






5.79 


Kr 


CELs 


4.21(-9)±1.24(-9) 


3.62±0.13 


+0.26±0.15 


+0.86±0.36 


+0.60±0.17 






>4.19 



References. — (1) Walsh et al. (1997) and (2) Dudziak et al. (2000) for 
& Liu (2005) for F, Sterling & Dinerstein (2008) for Kr, and Liu et al. (2004) 

*For Hen2-436, we adopted the CEL O abundance. 



TABLE 10 
P Abundances in 9 PNe 



Nebula 


[P/H] 


[Ar/H] 


[P/Ar] 


[C/Ar] 


Ref. 


1C5117 


+0.14 


-0.34 


+0.49 


+0.64 


(1),(2) 


NGC40 


-0.08 


-0.61 


+0.53 


+1.06t 


(3) 


NGC2392 


-0.65 


-0.21 


-0.44 


+0.34 


(4) 


NGC6210 


-0.23 


-0.19 


-0.04 


-0.12 


(5) 


NGC6826 -0.23 


-0.40 


+0.17 


+0.74 


(6) 


NGC6572 +0.34 


-0.27 


+0.62 


+0.21 


(7) 


NGC6741 


+0.58 


-0.01 


+0.59 


+0.52 


(8) 


NGC7027 


+0.71 


-0.25 


+0.96 


+0.98 


(9) 


Hen2-436 


+0.26 


-0.60 


+0.86 


+1.29 


(10) 



References. — (1) Rudy et al. (1991) for the P abundance; (2) Hyung et 
al. (2001) for the other elements; (3) see references in Table 9.; (4) Pottasch 
et al. (2008); (5) Pottasch et al. (2009); (6) Surendiranath & Pottasch (2008); 
(7) Hyung et al. (1994); (8) Hyung & Aller (1997); (9) Otsuka in prep.; (10) 
This work. 
* We adopted the CEL C abundance. 



Nordsieck 1977) with aniin=0.001 iim. and amax=0.25 (xm for 
amorphous carbon. For PAHs we adopted an oT^ size distribu- 



elemental abundances of Hen2-436. (3) Pottasch et al. (2003) for P and Fe, Zhang 
for the others of NGC40. 



tion with amin=0.00043 jjim and amax=0.001 1 jjm. We adopted 
an/?~^ dust density distribution. The dust size distribution and 
dust (and hydrogen) density distribution law are the same as 
those used in the Sgr dwarf galaxy PN BoBnl (Otsuka et al. 
2010). In BoBnl, Otsuka et al. (2010) used an am C + PAH 
grain model and estimated a dust mass of 5.78(-6) Mq and 
temperature of 80-180 K. 

In Table 1 1, we compare the predicted relative line intensi- 
ties with observed values, where /(H/3) is 100. The log g=6.0 
models give a better fit to the observations. Columns 3, 4, and 
5 list the observed values and the values predicted by mod- 
els A and B with log g=6.Q, respectively. Column 6 lists the 
observed values by Walsh et al. (1997) as a reference. In 
general, the predicted line and wide bands fluxes (except for 
[Oll]A3726/29 and 2MASS Ks) agree with the observations 
within 30% for models A and B. The large discrepancy be- 
tween the observed [O ii]A3726/29 line fluxes and the model 
might be due to the flux calibration uncertainty around 3700 A 
and the adopted monotonically decreasing density profile. 
Since the models did not include a density jump around the 
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Wavelength (/xm) 

Fig. 5. — The predicted SED from the P-I modehng (blue lines) by 
Cloudy. The observed data are indicated by the black lines or circles. In 
the inner small boxes we plot the radio data and the predicted SED. (upper 
panel) The predicted SED by CLOUDY considering amorphous carbon only. 
{lower panel) The predicted SED considering amorphous carbon and PAH 
grains (see text). 

ionization front which sometimes weakens the line intensity 
by collisional de-excitation, the models would predict larger 
[Oll]A3726/29 fluxes than the observations. The large dis- 
crepancy between the 2MASS Ks and the models might be 
due to the contribution from the H2 line. 

In Table 12, we list the derived parameters of the PN cen- 
tral star, ionized gas nebula, and dust. In Figure 5 we present 
the predicted SED from CLOUDY (blue lines). The upper 
and lower panels are the results of models A and B, re- 
spectively. In the inner small boxes we plot the radio flux 
densities (circles) from Dudziak et al. (2000) and the pre- 
dicted SED. Dudziak et al. (2000) measured flux densities of 
0.62/3.90/4.90 mJy at 1.46, 4.89, 8.40 GHz respectively. Al- 
though the predicted SED well matches the data in optical to 
radio wavelengths, it cannot explain the ^30 /im bump and 
the IRAS 100 /im data. The 30 /im bump might be from Mag- 
nesium Sulfide (MgS) grains, which are sometimes observed 
in C-rich PNe. At the present, since there are no available 
optical constants of MgS at UV wavelengths, we did not con- 
sider MgS in the SED models. Most of the IRAS 100 /im flux 
would be from the surrounding cold ISM because the spatial 
resolving power is ^2' at this band. Accordingly, our models 
failed to explain the flux in these bands. 

For Hen2-436, this work is the first successful estimate of 
the dust mass of 2.9(^) Mq in model A and 4.0(^) Mq in 
model B and temperatures of 100-150 K in model A and 100- 
200 K in model B. The derived dust-to-gas mass ratios in both 
models are slightly larger than the typical value in PNe hav- 
ing nebular- radius of 10'^ cm (--10"''; Pottasch 1984). The 
value is rather close to the LMC ISM (4x10"-'; Meixner et 
al. 2010). Lagadec et al. (2009, 2010) assumed a dust-to-gas 
mass ratio of 5 x 10"^ in radiative transfer modelings of C-rich 
stars in the Sgr dwarf galaxy. 

We feel that model B can better explain the physical con- 
ditions of gas and dust in Hen2-436, because the line fluxes 
predicted by model B give a better fit to the observations and 



TABLE 1 1 

The predicted relative line fluxes by the Cloudy models. 



Ion 


A 


/(A)„i„ 


' 'cloudy 


' V'^/cloiidy 


Walsh et al. 




(AZ/im) 




Model A 


Model B 


(1997) 


[0 II] 


3727 


9.25 


27.80 


25.60 


10.86 


[Ne III] 


3869 


53.57 


61.64 


56.87 


57.70 


[S II] 


4068 


2.07 


2.95 


3.12 


2.78 


IS III 


4076 


0.81 


0.95 


1.00 


0.98 


C II 


4267 


1.10 


1.19 


1.31 


0.97 


II 


4294 


0.04 


0.02 


0.03 




[0 III] 


4363 


13.90 


12.37 


10.84 


14.31 


He I 


4411 


5.91 


6.99 


6.53 


6.11 


lAr IVl 

L^-vi 1 V J 


4740 


0.73 


0.99 


0.82 


0.74 


[0 III] 


4931 


0.15 


0.11 


0.11 




10 ml 


4959 


265.07 


259.39 


262.20 


271.64 


10 ml 


5007 


792.07 


780.77 


789.22 


790.83 


[Ar III] 


5192 


0.06 


0.11 


0.13 




[Fe III] 


5271 


0.08 


0.24 


0.15 




[CI III] 


5518 


0.06 


0.04 


0.05 


0.21 


[CI III] 


5538 


0.16 


0.15 


0.17 


0.13 


[N II] 


5755 


0.72 


0.70 


0.73 


1.10 


He I 


5876 


19.56 


22.36 


20.73 


18.21 


[S III] 


6312 


1.63 


1.06 


1.14 


2.47 


[N II] 


6548 


3.59 


3.59 


3.80 


3.54 


[N II] 


6583 


10.92 


10.61 


11.22 


11.43 


He I 


6678 


4.92 


5.20 


4.88 


4.77 


[S II] 


6716 


0.49 


0.69 


0.66 


0.49 


IS III 


6731 


1.10 


1.42 


1.41 


1.10 


[Ar III] 


7135 


7.66 


8.03 


11.16 


6.48 


[Ar IV] 


7171 


0.04 


0.04 


0.03 




[Ar IV] 


7263 


0.04 


0.03 


0.02 




[O II] 


7323 


8.29 


16.84 


15.74 


6.24 


[O II] 


7332 


6.75 


13.59 


12.69 


7.03 


[Ar III] 


7751 


1.84 


1.94 


2.69 


1.22 


[CI IV] 


8047 


0.17 


0.21 


0.25 




C III 


8197 


0.05 


0.13 


0.10 




He I 


1.28 


1.13 


0.98 


0.93 




Hi 


1.28 


16.18 


15.63 


15.61 




[Ne III] 


15.55 


31.41 


29.87 


30.92 




2MASS J 


1.24 


51.90 


49.06 


50.57 




2MASS H 


1.66 


40.00 


27.37 


32.20 




2MASS Ks 


2.16 


98.50 


22.21 


23.81 




IRSB 


20.0 


2182.05 


2199.71 


2044.33 




IRSC 


30.0 


1004.30 


979.57 


922.92 





Teff is not so high as to produce strong emission from He II, 
which would originate in the stellar wind of Hen2-436. The 
predicted Teff and L^, are comparable to the values by Dudziak 
et al. (2000). At this time, there is no UV data of Hen2-436. 
The flux density profile in the UV region and the mid IR flux 
profile depend upon the central star temperature and luminos- 
ity as well as the dust composition. The predicted SEDs in 
the optical or longer wavelength regions by Models A and B 
do not show large differences, however a slight difference is 
found around 0.1-0.2 fj,m as can be seen in Figure 5. This 
discrepancy could be due to dust composition. Certainly in 
PNe, the bulk of the heating of the grains will be done by UV 
photons around 0. 1 /im, where the dust absorption coefficient 
peaks. To further enlarge our understanding, UV and also 3- 
10 /im spectra would be necessary. These data would enable 
us to estimate the C abundance using C [ill] AA1906/09 lines 
and the CEL C/0 ratio and also check for the existence of 
PAHs. 

In general, the elemental abundances of PNe have been es- 
timated based on gas emission lines only. Here, we can esti- 
mate gas and grain C mass through SED modeling. We can 



12 



Otsuka et al. 



TABLE 12 

The derived properties of the PN central star, ionized 
nebula, and dust by the p-i model. 



parameters 



log,L,/Z,0 
log,r,y/(K) 
loggCcm^ s"') 
composition 



parameters 



composition 



Central star 



Model A 



3.60 
5.07 
6.0 

[Z,F]=0, [Z] = 



-1 



Nebula 



Model A 



Model B 



3.57 
4.93 
6.0 

[X,F]=0, [Z] = 



Model B 



He: 11.15,C:9. 17,N:6.93, He: 11.08,C:9. 19,N:7.03, 



O 



0:8.29,Ne:7.43,S:6.15, 0:8.33,Ne:7.45,S:6.25, 
Cl:4.33,Ar:5.74,Fe:5.72 Cl:4.43,Ar:5.89,Fe:5.57 
others:[X/H]=-0.6 others:[X/H]=-0.6 



geometry 


spherical 


spherical 


RiJRout (") 


0.19/0.29 


0.18/0.26 


log.A'H(R,„) 


(cm-3) 4.75 


4.78 


logF(H,3) 


-12.02 


-12.02 


Mgas (Mq) 


0.05 


0.07 


dust comp. 


am. C 


am. C + PAHs 


Tdust (K) 


100-150 


100-200 


^dust (A^©) 


2.9M) 


4.0M) 


Mdust/A^gas 


5.92(-3) 


5.58(-3) 




TABLE 13 


Gas AND GRAIN CARBON 


.MASS IN HLN2-436. 




parameters 


Model A Model B 




gasC(M0) 


1.3(-3) 8.6(^) 




grain C(Mq) 


2.9i-A) 4.0(-A) 




gas + grain C(M0) 


1.6(-3) 1.3(-3) 




Af(gas + grain C)/Af(H)t 


1.8(-3) 2.3(-3) 



number density 



then investigate how much the C abundance increase if we 
include the contribution of grains. In Table 13, we present 
the gas and grain phase C mass in both models. It is remark- 
able that in model B, which is the best model at the present, 
^50 % of the gas C mass exists as grains. If the C abundance 
is measured using the sum of gas and grain C masses, the C 
abundance (logj,, A'(C)/A^(H)+12) would be 9.36 dex, which is 
~0.17 dex larger than the case of only gas C (9.19 in the case 
of gas C only, see Table 12). The grain C abundance might be 
considerable in PNe. 

4.2. Evolutionary Status 

In Figure 6 we plot the location of Hen2-436 and two Sgr 
dwarf galaxy PNe StWr2-21 and Wrayl6-423 (Zijlstra et al. 
2006) and the post-AGB He-burning evolutionary tracks with 
Z = 0.008 by Vassihadis & Wood (1994). Since aU these three 
Sgr dwarf galaxy PN have WR type central stars (Zijlstra et 
al. 2006), we assume that the central stars are hydrogen-poor. 
For Hen2-436, we find that predicted by models A and B 
has and uncertainty of 1 200 Lq , which is estimated from the 
H/3 flux and distance determination errors. The uncertainty 
of Tefc is ~ 10 000 K. These evolutionary tracks suggest that 
the progenitors are 1 .5-2 Mq stars, which end their lives as 
white dwarfs with a core mass of ^0.63-0.67 Mq. Among the 
above three PNe, Hen2-436 seems to be the youngest PN. The 
age of Hen2-436 is estimated from the evolutionary tracks to 
be ~3000 yr after leaving the AGB phase. Gesicki & Zijl- 
stra (2000) estimated the expansion velocity to be 14 km s"'. 
When adopting the value Rom = 0.26" (see Table 12) and the 
distance of 24.8 kpc, the dynamical age is estimated to be 
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Fig. 6. — The evolutionary tracks for He-burning stars with Z = 0.008 and 
the positions of three Sgr dwarf galaxy PNe. 



TABLE 14 

The OBSERVED CHEMICAL ABUNDANCES AND THE THEORETICAL 
MODEL. 



Nebula 


He 


c 


N 


O 


F 


Ne 


P 


S 


Ref 


Hen2-436 


11.09 


9.08 


7.38 


8.35 


5.69 


7.62 


5.72 


6.83 


(1) 


StWr2-21 


11.00 


9.00 


7.88 


8.53 




7.54 




6.93 


(2) 


Wrayl6-423 


11.03 


8.86 


7.68 


8.33 




7.55 




6.67 


(2) 


Model 


2.25 iW .; (Z = O.OOS) 


11.01 


9.07 


8.00 


8.54 


5.00 


7.95 


5.18 


6.85 


(3) 



References. — (1) This work; (2) Zijlstra et al. (2006); (3) Karakas 
(2010) 

~2200 yr, which is consistent with the evolutionary age. 

4.3. Dust Mass-Loss Rate 

If most of the observed dust is formed during the last two 
thermal pulses (^10000 yr, Vassiliadis & Wood 1993 for the 
case of a 2.0 M© initial mass and Z = 0.008), the dust mass- 
loss rate (Mdust) is estimated to be < 3 . 1 (-8) Mq yr~^ . Lagadec 
et al. (2010) estimated a Mdust of 0.97-2.52(-8) Mq yr"' 
for IRAS 18436-2849 in the Sgr dwarf galaxy. IRAS 18436- 
2849 shows similar metallicity and luminosity to Hen2-436. 
For a C-rich metal-poor star ([Fe/H] — 1) IRAS 12560H- 1656 
in the Sgr dwarf galaxy stream, Groenewegen et al. (1997) 
estimated Mdust to be 1.9(-9) Mq yr"'. Our estimated value 
is comparable to their values. Adopting a dust-to-gas mass 
ratio of 5.58(-3), the mass-loss rate M during the last two 
thermal pulses is estimated to be <5.5(-6) Mq yr"', which is 
very close to the predicted mass-loss rate during the last two 
thermal pulses (Vassihadis & Wood 1993 for the same case 
above). 

4.4. Comparison between the Observed Chemical 
Abundances and the Theoretical Model 

In Table 14, we present the observed elemental abun- 
dances of Hen2-436, StWr2-21, and Wrayl6-423 (Zijlstra et 
al. 2006) and predicted values from the theoretical models 
of Karakas (2010) for 2.25 Mq stars with initial Z=0.008. 
The abundances of Sgr dwarf galaxy PNe are measured us- 
ing gas emission lines only and their C abundances are de- 
rived from the ORL C lines, abundances from the models are 
the values at the end of the AGB. Karakas (2010) predicted 
that such stars end as white dwarfs with a core mass of 0.652 



Dust and chemical abundances of the Sgr dwarf galaxy PN Hen2-436 



13 



Mq, which is comparable to the core mass of Hen2-436 es- 
timated from the evolutionary tracks. In the model for 2.25 
Mq Stars with Z=0.008, a partial mixing zone, which pro- 
duces a '^C pocket during the interpulse period and releases 
free neutrons, is not included and the mass-loss law during the 
AGB phase by Vassiliadis & Wood (1993) is adopted. For the 
models with Z=0.008, Karakas (2010) adopted LMC compo- 
sitions from Russell & Dopita (1992). The accuracy of the 
predicted abundances by the models is on the order of 0.3 dex 
(A. Karakas in private communication). For the observed N 
and O abundances of Hen2-436, we adopted the values de- 
rived from the CELs. 

In general the observed chemical abundances of the above 
three PNe agree well with the model of Karakas (2010) 
within estimation errors, which suggest that all three of these 
PNe have evolved from w2.25 Mq single stars with initial 
LMC chemical compositions. The three Sgr dwarf galaxy 
PNe would have experienced evolution similar to LMC PNe, 
which is a remarkable finding. Some of the ''^N in Hen2-436 
might have been converted into ^^Ne by the double a particle 
capture process. The overabundances of '^F and ^^P in Hen2- 
436 imply that an extensive partial mixing zone was formed 
and the extra neutrons were released in the He-rich intershell. 

5. CONCLUSION 

We estimated elemental abundances in the Sgr dwarf galaxy 
PN Hen2-436 based on the archived ESO/VLT FORS2 and 
Spitzer/IRS spectra. We detected candidates of [Fll]A4790, 
[Krlll]A6826, and [Pll]A7875 for the first time, which indi- 
cates that these elements are largely enhanced. We found a co- 
relation between C and F, P, and Kr abundances among PNe 
and C-rich stars. The detections of F, P and Kr in Hen2-436 
support that F, P, and Kr together with C are certainly synthe- 
sized in the same layer and brought to the stellar surface by the 
third dredge-up. We detected some N ii and O ii ORLs and 
derived the ionic abundances from these lines. The discrep- 
ancy between O ORL and CEL abundances is > 1 dex. We 
constructed a SED model considering dust and estimated the 
initial mass of the progenitor to be ~ 1 .5-2.0 Mq with Z=0.008 



and the age to be ^3000 yr after the AGB phase. Hen2-436 
shares its evolutionary status with the Sgr dwarf galaxy PNe 
StWr2-21 and Wray 16-423. The observed elemental abun- 
dances of the three Sgr dwarf galaxy PNe could be explained 
by a theoretical nucleosynthesis model with initial mass 2.25 
Mq, Z=0.008, and LMC chemical abundances. The SED 
model predicted that >2.9(-4) Mq of carbon dust co-exists 
in the ionized nebula. Hen2-436 seems to have experienced 
evolution similar to LMC PNe. Based on the assumption that 
most of the observed dust is formed during the last two ther- 
mal pulses and the dust-to-gas mass ratio is 5.58(-3), the dust 
mass-loss rate and the total mass-loss rate are <3.1(-8) Mr.; 
yr"' and <5.5(-6) Mq yr"', respectively. Our estimated dust 
mass-loss rate is comparable to a similar metalUcity and lu- 
minosity Sgr dwarf galaxy AGB star. 
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